High temperature thermodynamics of fermionic alkaline earth atoms in optical lattices 
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We calculate experimentally relevant properties of trapped fermionic alkaline earth atoms in an 
optical lattice, modeled by the SU (N) Hubbard model. Our calculation is accurate when the tem- 
perature is much larger than the tunneling rate, similar to current regimes in ultracold atom exper- 
iments. In addition to exploring the Mott insulator-metal crossover, we calculate final temperatures 
achieved by the standard experimental protocol of adiabatically ramping from a non-interacting gas, 
as a function of initial gas temperature and final state lattice parameters. Of particular experimen- 
tal interest, we find that increasing N gives substantially colder Mott insulators, up to more than a 
factor of five for relevant parameters. This cooling happens for all N, fixing the initial entropy, or 
for all N < 20 (the exact value depends on dimensionality), fixing the initial temperature. 



Introduction. — The recent achievement of Fermi de- 
generacy and Bosc-Einstcin condensation in ultracold 
alkaline earth atoms [l[ opens up great opportunities 
in quantum information processing quantum sim- 
ulations Q, atomic clocks experiments and other 
precision measurements 0]. One fundamental prop- 
erty of fermionic alkaline earth atoms is their intrinsic 
S\J(N = 21 + l)-symmetry in the nuclear spin (I) de- 
grees of freedom, whilcbosonic alkaline earth isotopes in 
contrast have even-even nuclei and necessarily 7 = 00. 
Fermionic alkaline earth atoms loaded in an optical lat- 
tice are described by the SU(AT) Hubbard model where 
N can be varied for a single isotope from 2 to 21 + 1 < 10 
by selectively populating hyperfine levels. Cold atoms 
realizations of this model open up a range of exciting 
and exotic physics relevant to condensed matter: it is a 
simple limit of multiorbital models describing transition 
metal oxides, is important in theoretical generalizations 
of the Fermi-Hubbard model, and displays (frequently 
exotic) phenomena such as possible antiferromagnetism, 
supcrconductivity,nematic order, valence bond, and spin 
liquid phases 0, 0, Q • As a first step towards reaching the 
low temperatures necessary to observe these states, here 
we study the SU(N) Hubbard models high temperature 
Mott-metal crossover. One particularly interesting find- 
ing is that increasing N can lead to more than a five-fold 
decrease in temperature compared to N = 2. 

We calculate density and entropy profiles of lattice al- 
kaline earth atoms to second order in t/T, the tunneling 
rate over the temperature (see Eq. ([1])). This calculation 
is accurate for T t, regardless of the on-site interaction 
U. This includes the "unquenched Mott insulator (MI)" 
regime t < C/ that has been realized in SU(2) alkali 
gases 0. For the SU{2) spin-1/2 case, sophisticated, nu- 
merically intensive algorithms have yielded scries to tenth 
order |10| . For properties in the experimental regime, 
T > t in three dimensions, the second order expansion 
agrees quantitatively (< 1% error) with high order ex- 
pansions and dynamical mean field theory [llj. We also 
calculate final temperatures achieved by standard experi- 



mental adiabatic ramping protocols. Excitingly, our cal- 
culations show that applying the same protocols as in 
SU(2) experiments will generate colder, less compress- 
ible MI states as N increases, up to N ~ 10, with the 
most favorable N depending on dimensionality. 

Alkaline earth atoms in deep optical lattices are well- 
described by the SU (N) Fermi-Hubbard model [H 

H = —t fa,ifa,j + fasfa'.ifa',ifa,i (1) 

(ij),a i,a,a' 

where f a j is a fermionic annihilation operator destroying 
a particle of flavor a at site j, satisfying anticommutation 
relations {f a ,i,fa>,j} = S a ,a'Sij, J2(ij) indicates a sum 
over nearest neighbors, and a and a' are flavor indices 
that run from 1 to N. It is useful to note that the inter- 
action term may be rewritten as (U/2) Ei fhiifhi — 1), 
defining m, = E Q n a,i and n Qjl = f a J a .i- 

Atomic limit: experimental observables and interpre- 
tation in T/U <C 1 and T/U ^> 1 limits. — First we give 
results in the atomic limit (t = 0), the zero'th order 
term in the high temperature series expansion in t/T. 
Throughout, we will present explicit analytic results for 
the free energy density J 7 , from which the other consid- 
ered observables can be obtained by differentiating: the 
average filling and entropy per site are (m) = —dJ-fdfj, 
and s = —dT/dT, respectively. 

For the homogeneous system, the grand canonical free 
energy per site for t — is 



Jo = -Tlogzo 



(2) 



with the on-site partition function z = 

Em=o c, m e ~" eo(m) > where £ oM = {U/2)m(m-l)-ian, 
j3 = 1/T, ks = h = 1 throughout, and the binomial 
coefficient. The average filling is mo = { m )o where 
the expectation (O) of an operator O that depends 
only on the site filling is defined as the on-site expec- 
tation (O) = ^ E m C ' C 'm e " /3eo(m) and thc entropy is 
,s = log 2 + (1/T) (e (m)) . 
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C£, +1 ), with 



It is illuminating to consider the obscrvablcs in the 
T <C U and T > U limits. In the high temperature T > 
{U, fj,} limit the on-site partition function is z (T U) = 

Em=o C ™ = 2N > the density is m(T > U) = N/2, and 
the entropy density is s(T ^> U) = N\og2. Interpreting 
this, pairs of fillings to and TV — m are equally likely, for 
all m, so that the average filling is N/2, and each flavor 
is equally likely to occupy or not occupy each site, giving 
log 2 entropy per flavor per site. In contrast, in the low 
temperature MI limit defined by t <C T <C U and fi ^ 
Urn for all m, one term dominates z$ so that the partition 
function is zq = CT^e _ ^ e °( m ), the density is to chosen to 
minimize eo(m), and the entropy density is s = log(C^). 
Again, this may be simply understood: the filling is fixed 
at to and the entropy per site for the filling to MI state in 
this regime is just the number of ways to choose m of N 
flavors to put on a site. Finally, when t -C T <C U , but 
fx = Urn' — the metal boundary separating the to' and 
to' + 1 Mi's — the density is m = (to' + 1)/(1 + 1/TV), 
and the entropy density is s = log (C^, 
interpretations analogous to the discussions above. 

Cold atomic systems are confined in trapping poten- 
tials. To include their effects, we employ a Thomas-Fermi 
or local density approximation (LDA) [l2j ]. We take the 
trapped system's properties at a point in space r to be 
those of the homogenous system at a chemical potential 
/z(r) = fiQ — V(r), with V(r) the trapping potential. This 
is accurate when the potential varies slowly compared to 
the state's characteristic lengths, which is well satisfied 
in the regimes of present interest. For simplicity, and as 
an accurate description of most traps, we approximate 
V(r) = mui 2 r 2 /2 where m is the particle mass an u is 
the external potential's trapping frequency. 

Figure [T] shows the density to and entropy s profiles, 
and the formation of Mott plateaus at low temperatures. 
It also shows the effects of the tunneling corrections 
calculated later. At high temperatures T/U > 0.2 all 
curves are smooth and lack visible signs of Mott plateaux 
(not shown). Fig. [TJs left and right insets show the TV- 
dependence near m = 1 and m = N/2, respectively. Al- 
though the density profiles shown are at temperatures 
below the regime of validity t/T > 1, one finds that 
the theory is inaccurate only near chemical potentials 
halfway between the Mott shells • Working deeper in 
the approximation's regime of validity (smaller T) would 
lead to a smaller effect, invisible to the eye. However, the 
qualitative effects are the same, only smaller. 

Adiabatic loading. — We provide a theoretical descrip- 
tion for the standard protocol used to realize both 
bosonic [IH and fermionic 0, UH Mi's. This is essential 
to understand and optimize the process. Remarkably, we 
find that for the relevant TV, the final achievable temper- 
ature substantially decreases with increasing N. 

The procedure used to create Mi's is to first create a 
degenerate, weakly interacting gas without a lattice. A 
lattice is then ramped up to its final value. Ideally, the 
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FIG. 1. Observables as a function of distance to trap center r 
for N = 4 at T/U = 1/15, fi/U = 3.0, and in two dimensions. 
Top: density (solid: atomic limit; dashed: t/U = 0.1). Al- 
though the value of the t/T is not deep in the regime where the 
high temperature expansion is valid, such a large value was 
chosen so the effects would be visible in the density profile. 
Insets: density versus fi for TV = 2, 5, and 10, showing TV and 
m-dependence. Left inset: zoom around m = 1 shell. Right 
inset: zoom around m = TV/2 with density and /j, shifted by 
m s and fi a so that the Mott shells nearly overlap on the plot. 
Bottom: entropy (solid: atomic limit; dashed: t/U = 0.04; 
horizontal dashed: T <C U deep Mott and deep metal limits 
discussed in text). 



ramp is slow enough that the process is adiabatic. The 
adiabatic limit a ppe ars to be approached in recent boson 
13l UM i and with less than ~ 50% increases 



experiments ^ 

of entropy in some SU(2) Fermi experiments [tllflij]. The 
limits of adiabaticity are ill-understood and beyond the 
scope of our present work. 

In the adiabatic limit entropy is conserved, and given 
the initial state's particle number and entropy one can 
determine the final temperature by matching the parti- 
cle number and entropy to the initial state. The initial 
state entropy may be controlled and measured only in- 
directly, however, through the temperature. Thus, we 
must first determine S of the initial state from the initial 
temperature T±. 

For any initial state sufficiently cold to reach a MI, 
the initial gas will be deeply degenerate, T <^ \i. For 
large numbers of particles, the harmonic trap can be 
treated as having a continuous density of states v(e) — 
(e/o;) d - 1 /(w(d - l)!)©(e) with 6(e) the Heaviside step 
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function. Using this, the total particle number and 
entropy of a d-dimensional trapped system are M — 
f (£)" and S = (a)"" 1 to lowest order in Tj \i. 

Note that to the same accuracy as this approximation, 
the Fermi temperature is Tf = \i- 

Figure [2] overviews the results of our adiabatic loading 
calculation for a two dimensional system. Although spe- 
cific values depend on microscopic parameters such as the 
scattering length a s , lattice spacing a, and trap length 
£, the qualitative findings discussed below are indepen- 
dent of these. Fig. [2] shows the final temperature T/U 
after adiabatic loading as a function of the initial tem- 
perature Ti for several TV, using experimentally relevant 
parameters at fixed final lattice depth. Fig. [2] (inset) il- 
lustrates the effect on density profiles of cooling obtained 
by increasing TV. For reference, the coldest alkaline earth 
gases in the weakly interacting regime are at rou ghly 
T/fi = 0.14 for 173 Yb (TV = 6) and T/n = 0.26 [3 
for 87 Sr (TV = 10) 0. Also, we have found that in the 
deep lattice regime, the final rescaled temperature T/U 
is independent of lattice depth Vq within a few percent. 

One of our most remarkable findings is the final tem- 
perature's TV-dependence: increasing TV produces colder 
MPs with initial parameters fixed in an experimentally 
realistic way. In particular, this result holds when the ini- 
tial temperature of the weakly interacting gas is fixed, as 
would hold if the gas were sympathetically cooled with 
another species to prepare the initial state — an even 
more favorable situation is considered below. A corol- 
lary is if adiabaticity can be maintained, since SU (TV) 
systems can be produced with initial temperatures com- 
parable to current SU (2) systems MPs with T <C U and 
even T ~ t are well within reach for ST/ (TV) fermion ex- 
periments. 

To understand our findings, we consider the TV- 
dependence of S in the initial and final states. In the 
initial state, S at fixed M scales as S{ oc N x / d . For the 
MI at T/U <C 1, the entropy scales as logC/^. Consider- 
ing the m = 1 MI shell for simplicity, Sf oc log TV. Thus, 
naively one may expect the initial state's entropy to grow 
faster with TV than the final state's at fixed initial and 
final temperatures, but this is true only for very large 
TV: in practice, for d = 3, the power law grows slowly 
compared to the logarithm, and for TV < 20 the ratio 
of MI entropy to metal entropy increases with TV. This 
results in cooler states with increasing TV in this regime, 
as observed in Figs. [5] Similar effects persist in lower di- 
mension, but to smaller TV (roughly 7 and 3 in d = 2 and 
1, respectively). For m that scale with TV, e.g. m = TV/2, 
the situation is even more favorable: one finds that the 
MI entropy is proportional to TV and thus for any d > 1 
the final states get colder with increasing TV for all TV. 
A similar argument explains the reversal of this effect at 
higher temperatures (see Supplementary information). 

While we observe this dramatic cooling when fixing 



T, as appropriate, say, for sympathetic cooling, the ac- 
tual situation may be even more favorable. Other cooling 
procedures may fix the initial gas's entropy or T/fx inde- 
pendent of TV. In this case the MI will get colder with 
increasing TV for all TV, even for m = 1. A simple candi- 
date method is direct evaporative cooling of the fermions, 
with an added benefit that Pauli blocking becomes less 
important as TV increases. 




Ti/co 

FIG. 2. Left: Adiabatic loading: final temperature T/U as 
a function of initial temperature T/uj for TV = 2,4,6 (top 
to bottom at low temperatures) at lattice depth Vo = WEr, 
with En = fi 2 Tv 2 /(2ma 2 ) the atom recoil energy, in two di- 
mensions. Parameters are fixed to experimental values for 
173 Yb[jJ]: we find U/uj = 50, a/ 1 = 0.4, and M = 5 x 10 4 
using standard optical lattice results for U [l7j|- The experi- 
mental initial gas temperature is T/ui = 5.2 (T/fi = 0.14 for 
N = 6). Parameters are roughly the same for other alkaline 
earth atom experiments. Right: Density profiles after adia- 
batic loading for T/u = 7.3 (T/a = 0.14,0.17, and 0.20 for 
TV = 2,4,6, respectively). 



Non-zero tunneling. — Here we incorporate tunneling 
by a second order expansion in t/T, equivalent to a fi- 
nite temperature t/U expansion. Analogous results were 
recently applied to SU(2) Fermi gases [ll|, where this 
second order series is quantitatively accurate (< 1% er- 
ror) when T > t for fillings in the metallic region (e.g. 
77i ~ 1.3) and for much lower t emp eratures at other fill- 
ings (e.g., m ~ 0, 1, and 2) [l8l.ll9|. Present experiments 
are in this regime, although SU (2) experiments have be- 
gun to reach temperatures T <t. Longer series maintain 
agreement for T > t, but, being divergent, scries of all 
lengths give unphysical results for T <t. Longer series 
may nevertheless be useful in understanding low temper- 
ature physics by analyzing their analytic structure |10j . 

We find the free energy density to 0[(t/T) 2 ] is 
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where we define 
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Eq. ©'s 



with gr m4 _ w m ^ z ^ ( m _ z+1 )2 e ^(« (m)+ eo (i)) 
bracketed term comes from Pauli blocking. 

Figure Q] shows tunneling's consequences for a two di- 
mensional system. In the density profiles, we see that it 
reduces the Mi's size, as expected. In the entropy pro- 
files, tunneling slightly increases the MI entropy while 
significantly decreasing the metal entropy: in the MI, 
tunneling reduces the excitation gap, increasing the en- 
tropy, while in the metal, the increasing bandwidth low- 
ers the low energy density of states, decreasing the en- 
tropy. We also observe (Fig. Q] insets) that although 
the magnitude of the tunneling corrections depends only 
weakly on m, the iV-dependence of the tunneling correc- 
tions depends more strongly. For m = 1, there is signif- 
icant iV-dependence of the tunneling corrections, while 
there is little in the m = N/2 shell. 

We finally have considered the effect of tunneling on 
the adiabatic loading procedure, and find that tunneling 
increases the final temperature. This may be understood 
by considering Fig. [I] tunneling reduces the total entropy 
in the trap as a consequence of a small increase in MI 
entropy and a relatively larger decrease in metal entropy. 
The effects are small when t/T is small, the limit in which 
our calculations are valid. For example, an N = 5, d = 
2 system with U/uj = 50 and a/£ — 0.13 initially at 
T j fx = 0.2 adiabatically loads to a MI with temperature 
T/U = 0.12 for t/U = and T/U = 0.14 for t/U = 
0.1. Tunneling corrections somewhat decrease the cooling 
effect of increasing N. 

Conclusions and discussion. — We studied the mctal- 
Mott insulator crossover using a high-temperature series 
expansion technique, up to second order in t/T. We cal- 
culated the density and entropy in the atomic limit and 
quantified how tunneling reduces the size of the Mott 
insulating region, its flavor number (iV) and filling (m) 
dependence, and how it increases the Mott entropy while 
decreasing the metal entropy. 

Additionally, we studied the standard experimental 
protocol used to realize Mott insulators and showed that 
the final temperatures significantly decrease with increas- 
ing TV. Our explicit calculations were in two dimensions, 
but similar results hold in other dimensions, becoming 
weaker in one dimension and more dramatic in three di- 
mensions. This is an encouraging result for attempts to 
reach large-iV Mott insulators, which paves the way to 
study exotic physics. The iV-dependence of low temper- 
ature T <C t behavior — e.g., antiferromagnets or other 
ordered or exotic states — remains to be investigated. 



Three body losses might limit the observability of the 
shell structure with m > 2. Nevertheless, based on 
Ref. [H's theory we find, at least for 87 Sr, that even fill- 
ings m < 5 may live sufficiently long to explore their 
many-body physics (see Supplementary material). 
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OTHER OBSERVABLES 

In the main text, we concentrated on the density and 
entropy. This was largely because the density is the most 
common observable, while the entropy controls the tem- 
perature attained during adiabatic loading. However, 
other interesting quantities may be measured. 
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FIG. 3. Top: compressibility function of distance to 

trap center r for N = 4 at T/U = 1/15, fi/U = 3.0 (solid: 
atomic limit; dashed: t/U = 0.1). These parameters corre- 
spond to the density profile shown in the main text. The 
dips in the compressibility at the metal center are artifacts 
of applying our approximation for a t/U that is larger than 
where the approximation remains valid. This was to ensure 
that the effect of tunneling would be visibile in the density 
profile, in the main text. Bottom: atom number distribution 
showing Pi, P2, P3, and P4 for the same parameters as the 
compressibility, in the atomic limit. 



Here, we discuss the compressibility k = dn/dj-i = 
—d 2 J-/d[i 2 and the on-site atom number distribution P m , 
the probability of m particles occuring on a site. The 
former is simply a derivative of the density (or the free 
energy) and is thus readily available from our calculation. 
The latter is not simply a derivative of the free energy 
we compute, without adding extra terms in the SU(N) 
Hamiltonian, and so we only consider it in the atomic 
limit. In this limit, one obtains P m = (l/z )C ? ^e"^ eo( ™ ) 
using the same notation as in the main text: eo(m) = 
(U/2)m{m - 1) - ixm and z Q = £™=o C^e~^ m \ It 
is worth noting that once one knows P m , one knows all 
moments of the on-site density: (m), ("i 2 ), . . . However, 
note that although for the SU (2) case the double occu- 
pancy is uniquely determined once one knows the density 
and number fluctuations (n 2 ) — (n) 2 , the same is not true 
for the SU(N) gas. 

Figure [3] shows k and P m as a function of the distance 
to trap center for the same parameters as Fig. 1 (bot- 
tom) of the main text. For these parameters, the com- 
pressibility is nearly zero in the Mott regions and peaks 
in the metal as expected. Tunneling reduces the magni- 
tudes of the peaks and slightly increases the magnitude 
of the compressibility in the Mott insulator as expected. 
At larger temperatures T > 0.3C7, this structure in the 
compressibility becomes quite washed out. The atomic 
limit number distribution shows that in Mott regions, a 
single filling dominates the number distribution, while in 
the metal region there are two relevant fillings whose fill- 
ings change with r to smoothly connect the Mott regions. 
At higher temperatures T > Q.3U, the Mott regions will 
have significant probability to occupy fillings other than 
the Mott filling and the metallic region will have signif- 
icant probability to occupy fillings other than the two 
fillings that are relevant at low temperature. At suffi- 
ciently high temperature, all of the structure is washed 
out. 



WHY, AT HIGH TEMPERATURES, HEATS 
RATHER THAN COOLS WITH INCREASING N 

Figure 2 of the main text shows that at low temper- 
ature, the effect of increasing N for experimentally rele- 
vant N while fixing the initial temperature is to generate 
colder clouds. However, it also shows that the trend is 
reversed at higher temperatures. This may be somewhat 
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counterintuitive, as at very high temperatures the en- 
tropy should become very high even in the deep lattice. 
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FIG. 4. Top: density profiles for N = 2, 4, 6 (bottom curves 
are N = 2 and N = 6, colored blue and gold, respectively, 
and are extremely similar; top curve is N = 4.) of final states 
after adiabatic loading for Tt/uj — 14.7 and other param- 
eters corresponding to Fig. 2 of the main text. Note the 
absence of a Mott regime at these high final temperatures. 
Bottom: entropy profiles for N = 2,4,6 (top/blue, mid- 
dle/dashed/purple, bottom/gold respectively), rescaled by 
r 2 /JV 1/ ' 3 , of final states after adiabatic loading for Ti/cu = 14.7 
showing that the entropy scales as iV 1//3 at these high tem- 
peratures rather than the more favorable log N obtained at 
lower temperatures. (The r 2 is the surface area factor that 
multiplies the entropy density in the integral to obtain the 
total entropy, and is included simply so one may more easily 
see what regions of the system contribute most to the total 
entropy.) 

To understand this behavior, we first note that at the 
higher temperatures the density profiles are quite differ- 
ent than at low temperatures: they have lower central 
density and are more spread out. Fig|4] illustrates this 
effect. In particular, there is no central Mott regime in 
the final states and it turns out this reduces the entropy's 
dependence on JV. In particular, we recall that at low 
temperatures, the final state entropy grows with log N, 



while the initial state grew with TV 1 / 3 in three dimen- 
sions; for the N < 20 of interest, the former grew faster 
than the latter, leading to the cooling. In contrast, Fig.0] 
illustrates that with the lower density clouds produced 
at high temperature, there is no Mott region and the en- 
tropy always scales roughly as iV 1 / 3 rather than log AT, 
eliminating the cooling effect seen at low temperatures. 
The details of this heating will depend on the central 
filling, and will presumably usually occur at higher tem- 
peratures for higher central fillings. 

LOSS RATES 

As a final note we mention particle losses. The dom- 
inant loss source will be three-body losses with recom- 
bination rate 7 (m(m — l)(m — 2)). For filling m = 1 
and m = 2 Mott insulators, the lifetimes will be quite 
long, limited only by losses from number fluctuations to 
higher fillings for all alkaline earths. For larger fillings, 
we estimate loss rates based on the theory of Ref . [l[ , av- 
eraging interference of the two main recombination path- 
ways. This is accurate to a few percent for both small 
a and large a and overestimates losses in between. For 
the 87 Sr experiments, we find 7 = 0.03s -1 in a very deep 
Vq = 35-Er lattice, with slower loss rates in weaker lat- 
tices. Correspondingly, the m = 3 shell in the atomic 
limit (no number fluctuations) has a lifetime T3 ~ 6s. 
For 173 Yb experiments, due to the roughly doubled scat- 
tering length we find six times faster loss rates at the 
same lattice depth (in recoil units) Vb / Er and the same 
lattice spacing as in the previous Sr estimates. 

However, this comparison is complicated by three 
points. Firstly, the energy scales in the Hubbard Hamil- 
tonian also increase, by a factor of about two due to the 
doubled scattering length. Hence, the relevant quantity 
— the ratio of the loss rates to those energy scales — is 
only larger for 173 Yb by a factor of three fixing Vo/Er, 
working at the same lattice depth and spacing. 

Unfortunately, there is a second consideration which 
makes things much less favorable for current experi- 
ments. Current 1 Yb experiments are performed in lat- 
tices with wavelength shorter by a factor of about two, 
which increases the loss rates enormously (~ 2 6 ). This 
factor comes because, for fixed Vo/Er, halving the lattice 
spacing in three dimensions increases the density by 2 3 
and the on-site losses go like maximum density squared 
(an integral over space of the density cubed J drp 3 (r)). 
This is somewhat compensated by an increasing Hamil- 
tonian energy scales. The interaction increases propor- 
tional to the maximum density (~ J drp 2 (r)), by a factor 
of ~ 2 3 . This gives an overall loss rate, relative to rates 
governed by the Hamiltonian, a factor of 8 larger than 
the estimate above, at the same Vo/Er. This still may 
be tolerable; the m = 3 Mott insulator has a lifetime 
of about 40ms and in deep lattices the lower shells have 
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very long lifetimes. Furthermore, this constraint to use 
a much shorter lattice spacing is not fundamental. 

Moreover, a third consideration shows that these esti- 
mates are somewhat too pessimistic. We assumed that 
Vo/Er was fixed, but both the increase in scattering 
length a s and decrease in lattice spacing a decrease the 
Vo/Er required to obtain the same t/U, thus counterbal- 
ancing the increased loss rates somewhat. The details of 



the balance of these effects depend on the regime of the 
phase diagram being explored. 
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